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SUMMARY 


In the method of Helmholtz-Klrchhoff for separated flow past a flat 
plate (normal to the stream) the separation velocity and the "base pres- 
sure" are fixed at the free -stream values. In the present treatment a 
modification is introduced to allow arbitrary separation velocity and 
base pressure, so that values more in conformity with experiment may be 
chosen. The solution depends then on the single (base-pressure) param- 
eter k. When k is suitably chosen, the drag and the details of the 
potential flow near the plate agree well with experiment. The compu- 
tations depend on a particular choice of free-streamline hodograph, which 
has the feature that it gives a definite wake width for every value of k. 
In this way the wake width is correlated with the drag. 

The same ideas are applied to work out the free-streamline flows for 
a circular cylinder and a wedge of 90 ° vertex angle. 


INTRODUCTION 


It is remarkable that the problem of flow past bluff bodies, one of 
the earliest to receive attention, is not yet understood. The early 
Investigators already had very good insight into the problem; and, 
although considerable experimental and some theoretical Information on 
the matter has been collected since then, there has been little essen- 
tial progress toward a theoretical formulation. One of the methods used 
to attack the problem was that of the free-streamline theory, introduced 
by Helmholtz and extended by Kirchhoff and many others. Kirchhoff’s 
example of flow past a flat plate, normal to the stream, is well-known. 
The theory is based on the observation that for the configurations in 
question the flow separates from the body, leaving behind it a wake and 
creating a pressure drag quite distinct from that due to shearing forces 
on the surfaces. The main aim of the theory is to find the "free stream- 
lines" defining the wake, outside which the flow is potential, and to 
compute the resulting pressure drag. 

In the Kirchhoff theory there is a basic assumption which results 
in a considerable loss of reality. This is the assumption that the 
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velocity on the free streamline at separation is equal to the free-stream 
velocity U^. The pressure at the separation points, and on the "base" 
of the body, behind the separation points, is then equal to the free- 
stream pressure. This is not in agreement with experience, which shows 
that the base pressure is actually always lower than the free-stream 
value and that the drag is higher than that calculated by the Kirchhoff 
theory. However, the Helmholtz-Klrchhoff assumption is attractive, in 
that the velocity all along the free streamline may be considered con- 
stant (Uoq) up to infinity, which leads to a simple formulation in the 
hodograph plane. The theory has been applied to many shapes other than 
the normal flat plate; in fact, there is a rather expensive literature 
on the subject. However, in all cases (with the possible exception of 
cavity flows) the computed results fail to agree with experience, the 
discrepancy being mainly due to the basic assumption about the separation 
velocity. 

Clearly, if any progress is to be made with the free-streamline 
theory, it is necessary that the separation velocity be allowed to assume 
values different from Uoo. The modification required in the theory may 
be summed up in a single parameter k, which defines the separation 
velocity U 8 = kUoo and gives a base-pressure coefficient Cp 3 = 1 - k^. 
The base-pressure coefficient is always less than zero, corresponding to 
k > 1. Of course, it is not known what value of k should be assumed, 
but this is a problem that cannot be determined by the free -streamline 
theory. It must come from other considerations, principally of the 
dynamics of the wake. 

One might doubt that the free -streamline theory is applicable at 
all, particularly in view of the last remark. To sum up the evidence for 
and against it, the following experimental observations may be useful: 

(l) The discontinuity surfaces, or free streamlines, idealized in the 
theory, are well approximated by the actual shear layers that exist in a 
real fluid, for sane distance downstream of the separation points (ref. 1, 
p. 555)* (2) On the back of the obstacle, downstream of the separation 

points, the distribution of pressure is remarkably constant for almost any 
form of bluff body, even the extreme example of a flat plate inclined at 
small angle of attack (ref. 2 or ref. 3> P- 879)- This means that the 
velocity at the two separation points is the same, a fact which Is rather 
essential to the theory. (3) The shear layers do not continue far down- 
stream as assumed, but "roll up" to form vortices, alternately on each 
side. This vortex formation occurs behind all bluff bodies, provided 
there Is no interfering barrier between the separated shear layers, at a 
frequency which is characteristic for each body shape. Fage and Johansen 
(ref. 4) noted that the vortex frequency for different bluff bodies could 
be correlated by expressing it as a dimensionless frequency based not on 
the body dimension but on the distance d' between the shear - layers, 
measured at the section where they become parallel, before "rolling up." 
Whether or not the vortices are formed, the idea of free streamlines 
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extending to infinity is unrealistic, for the shear layers diffuse rapidly. 
Therefore, it appears rather fruitless to he too concerned with the details 
of the free streamline at infinity; the main aim should be to obtain a 
solution which is valid near the body. A more realistic way to formulate 
the problem is as follows: 

The flow past a bluff body is considered in two parts. Near the 
body it may be described by the free -streamline theory, provided that the 
parameter k is properly chosen. The description of the wake farther 
downstream must come from other considerations; if it can also be obtained 
in terms of the parameter k, then a complete solution may be found by 
joining the two parts. 

Even if it is not possible to complete the solution in this way, on 
theoretical grounds, the results of the free-streamline theory should be 
quite useful, particularly for correlating various bluff shapes experi- 
mentally. For instance, if from a study of one or two cases it can be 
determined how k depends on wake breadth, condition of the separated 
shear layers, and so forth, it may be possible to choose the appropriate 
value of k for any other bluff body. As another example, the relations 
between "bluffness," drag coefficient, wake breadth, shedding frequency, 
wake energy, and so forth, might be easily classified. 

The Kirehhoff solution for the normal flat plate is reviewed. Then 
it is shown how a more realistic solution may be obtained by allowing 
arbitrary base pressures. The examples of a wedge and circular cylinder 
are also worked out. 

The research was conducted at GALCIT under the sponsorship and with 
the financial assistance of the National Advisory Committee for Aeronautics. 
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streamwise distance to section where streamlines become 
parallel 


drag coefficient 
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C P 

Cps 

<5 


pressure coefficient 

value of Cp at separation point and on base 
average of Cp (over a surface) 



d 

d' 


S*G,H,I,J 



2k 


k - 1 S / U oo 

L 

n 

q. 

^S 


average value of Cp on front of plate or wedge 

breadth of any cylinder measured normal to stream 

distance between parallel free streamlines 

functions defined for convenience of calculation, 
appendix B 


length of wedge measured along a side 
wedge-angle parameter, 2a/it 
magnitude of velocity (dimensionless) 
value of q. at separation point 


R 

t 


u,v 

w 

X 


y 

z 


Reynolds number 
intermediate mapping function 

magnitude of free-stream velocity (dimensionless) 

components of velocity 

complex stream function, 9 + ii|f 

streamwise coordinate in real plane 

coordinate normal to flow 

dimensionless complex coordinate, x + iy 
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0s 


s 

0 

9 

v 

P 

T 

X 

* 

ft 

fl k 

ft 1 


half -angle of wedge 

angular position on cylinder circumference, measured from 
stagnation point 

angular distance to separation point on circ ular cy lin der 

1 1 ie 

inverse of complex velocity, — = ~ e 

direction of flow 
velocity potential 

ifl 

complex velocity, qe 
density 

Levi Civita plane 
intermediate mapping function 
stream function 

intermediate mapping function, log e £ 
intermediate mapping function for k ^ 1 

"shaping" term in ft 

KIRCHHOFF PROBLEM 


The problem made famous by Kirchhoff was that of the flow past a 
flat plate set normal to the stream. It will be useful to review his 
solution, following essentially the notation of Lamb (ref. 3, p. 99). 
In sketch 1 the z-plane is the plane of the actual flow. The solution 
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consists, as usual, of mapping z onto the jplane of the complex poten- 
tial w = cp + iijr. This is accomplished through the medium of the complex 

holograph or v-plane, where v=|^ = u-iv = qe“ i6 is the complex 

velocity. It is with this plane that the elegance of the TTe.1mho.lt 7 — 
Kirchhoff method is realized. The essential assumption is that every- 
where on the free streamlines SI and S' I the velocity is eq ual to 
the free-str.eam value IT,*,. Then in the holograph plane the free stream- 
line is simply the circle |v| = (=1 after normalization). Sec- 

ondly, the trace of the plate SS' in the holograph plane is known, 
since the flow direction there is constant . Thus the boundaries of the 
flow in the holograph plane are known a priori and are of simple form, 
so that the mapping to the w-plane is easily accomplished. 

In practice it is convenient to use £ instead of v, where 
£ = — = — Thus £ gives the true flow direction and the 

V 4 U.W 

reciprocal of the velocity magnitude at the corresponding point in the 
physical plane. Once the mappings are known, the solution is given com- 
pletely by 



q 2 = — = q 2 (w) 

Ul 2 


( 1 ) 


p 

where q determines the pressure, since Bernoulli's equation may be 
used to evaluate the pressure coefficient 


= 1 - 


Prom relations (l) the pressure everywhere in the flow field may be 
computed. 


In particular, at separation and all along the free stre aml ine 
Cp ■= Cp S - = 0, since q s = 1. The pressure coefficient in the wake and 

on the back side of the plate is also zero. The drag then Is simply due 
to the excess of pressure on the front; its value in the Kirchhoff example 
is C D = 0.88. 

Now actual experience shows that the drag is considerably larger, 
being of the order Cp = 2; the increase Is due mainly to suction on the 
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back of the plate, Cp g < 0. Corresponding to this, the velocity on the 

free streamline at separation is higher than the free- stream value. If 
qg = kLk>, k > 1, then Cp S = 1 - k2, and the drag contribution from the 

back side is k^ - 1. 

The difference between the computed and the actual drag is a serious 
discrepancy in the Kirchhoff solution, as well as in the many solutions 
which have been worked out for various body shapes using the same tech- 
nique. In other respects, the flow in the imme diate vicinity of a flat 
plate normal to the flow does resemble the conditions anticipated in the 
theory, as already stated in the "Introduction," so that the free-streamline 
theory need not be abandoned. It would appear necessary to modify the 
theory only to the extent of allowing the velocity at separation to assume 
an appropriate value qg = kU^. 

Such an adjustment may be made, in fact, by introducing the holograph 
of sketch 2. Here, the velocity at separation is allowed to be q s = kU^ 


© 

B* 

B 


Sketch 2 

and to remain at this value along the free streamline (circle in the 
holograph plane) until the latter becomes parallel to the free stream 
(point B in the holograph plane). At Infinity the flow must have returned 
to the point I, so that the free streamline in the holograph Is simply 
drawn by joining BI, giving the notch shown in sketch 2. The singularity 
(doublet) is still at I, as in the Kirchhoff hodograph, which is shown by 
a dotted line for comparison. 

This "notched hodograph" Is convenient, for it may be easily mapped 
onto the w-plane ; but it also approximates experiments, since the pressure 
on the free streamline does tend to remain constant for seme distance 
downstream of the separation point. The solution for flow past a normal 
flat plate, using this hodograph, is worked out in the next section. 
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NOTCHED HODOGRAPH 


Sketch 3 shows the planes needed to map the flow from the z- to the 



Sketch 3 



B B' 
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B' 
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t 
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w-plane , under the assumption that in the hodograph plane it is like 
sketch 2. To make the radius of the circle in the v- and £ -planes equal 
to unity, the velocity at infinity is set at Uoo = l/k; that is, q s = 1. 
The transformations are: 



which is the Joukowsky transformation, 



+ X2 
+ 1 


t 
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where 


h = 



and 

v = l/t 2 


Solving for £ gives 



where 


a 


= k 2 + 1 
k 2 - 1 



the constant of integration having been determined from, the condition 
z = 0 at w = 0. Relations ( 2 ) and (3) give the complete solution. 
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Of particular interest are points on the plate and along the free stream- 
line, which correspond to real values of w(=cp), so that the expressions 
there are somewhat simplified. 

To locate these positions in the physical plane z = x + iy, one 
has from equation (3): 




Equations (4), ( 5 ), and ( 6 ) locate positions on the plate, on the free 
streamline from S to B, and on the free streamline from B to <», 
respectively. At B the free streamlines become parallel, separated by 
the distance d' , which may he found by putting cp = a2 in equations ( 5 ) 
or (6): 




(T) 
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Also, the downstream distance to B is 


b , ( £±if + leg. 


/ k^ + lY 


6 k 2 - 1 


(Ta) 


The pressure coefficient may he evaluated from. 


k2 


C L « 1 - k 2 q 2 = 1 - 
P Ul2 


U 

The values of 1 1, | 2 are found from equation (3), which gives, for the 
positions corresponding to equations (4) , (5), and (6) above, the 
expressions 


( 8 ) 


6 - 




HI -1 




0 ^ cp < 1 


1 S 9 < a 2 > 


cp > a 2 


( 9 ) 


To find the drag coefficient, the average pressures on the front anrt 
back of the plate are computed: 


C =- i T / C dy=^r (l- k 2 q 2 )i 
P a J 0 p ^ dj 0 V * J q 


dep 


( 10 ) 


where use has been made of the relation 
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on the plate. Substituting from equations (9) gives for the average 
pressure on the front 


P — 1 

0 pf " a 


frit H S . tan-1 


2k 


k 2 -1 


- 2(k2 - 1) 


( 11 ) 


On the back the pressure is constant, so the average pressure there is 
simply 


°ps = 1 - k2 


( 12 ) 


Finally, the drag coefficient is 


°D = C pF “ Cps ( 13 ) 

The breadth of the plate d is easily found from equations ( 4 ) by setting 
cp = l: 


d = 



2k , k 2 + i , -1 k 2 - 1 

+ — tan — * 

k 2 + 1 k 2 - 1 a , 




The expressions derived above conrpletey describe the flow, for any value 
of the parameter k. For k = 1 they reduce to the Kirchhoff solution. 
In real flows, however, k > 1 and the problem remaining is to choose 
the correct value. This problem will be returned to later (see the sec- 
tion "Discussion"). First it is necessary to test the validity of the 
notched hodograph by comparing it with experimental results. For this, 
there are available same excellent measurements by Fage and Johansen 
(ref. 2). 


COMPARISON WITH EXPERIMENTS AND WITH ANOTHER HODOGRAPH 


In one of the cases studied by Fage and Johansen the velocity at the 
edge of the plate and the pressure on the back correspond to k = l.jh. 
Using this value in equations (8) and (9), "the pressure distribution on 
the front of the plate may be computed and compared with the measured 
values. Figure 1 shows that the agreement is excellent. 
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Now it is not at all clear that some other holograph might nob 
give eq ually satisfactory agreement, especially since the "end. points" 

of the pressure distribution ^Cp = 1 and Cp S = 1 - k^ are fixed once 

k has been chosen. Any other holograph curve joining the points S 
and I (sketch 2) will give the same base pressure coefficient Cp S 
but will otherwise change all the values that have been computed above. 
To investigate the sensitivity of such a change, SIS' was chosen to 
be a slnooth curve, as shown in sketch defined in such a way that its 



inverse (in the £ -plane) is an ellipse. This permits an easy trans- 
formation, the details of which are given in appendix A. The pressure 
distribution on the plate calculated for the "elliptical" hodograph is 
also shown in figure 1. It does not agree with the measurements so well 
as does the case of the notched hodograph. 

A much more sensitive comparison is given by the pressure distribu- 
tion along the free streamline, as shown In figure 2(a) for the two 
cases . The superiority of the notched hodograph Is Indicated by a com- 
parison with the experimental curve (measured for this purpose), which 
shows that the pressure tends to remain constant at first, as antici- 
pated. Figure 2(b) shows the streamlines computed for the two holographs 
and for the Kirchhoff case. The shaded region is the actual shear layer, 
measured by Fage and Johansen, which the free stre aml ines are intended to 
approximate. 

An unusual and very useful feature of the notched hodograph Is that 
the free streamlines become parallel at some section B-B 1 . In this way 
a definite value of the wake width d T is defined for every value of k. 
The definition of a wake width opens some new possibilities which will be 
taken up later. 
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Figure 3 shows how the drag coefficient Cp and the wake width 
d'/d depend on the base-pressure parameter kj the calculated values 
are given in table I. 

After the above method was worked out, the author 1 s attention was 
drawn to another method, which was given by Riabouchinsky (refs. 5 and 6 ) 
and which accomplishes essentially the same thing. Riabouchinsky intro- 
duces, further downstream, a second plate which is the mirror image of 
the first one. The two plates and the two streamlines which join their 
corresponding edges enclose a region, or cavity, which is assigned an 
arbitrary pressure. This is then the base pressure, as well as the con- 
stant pressure along the free streamline, and may again be specified by 
the parameter k. For a given k the outer, potential flow and the shape 
of the cavity are given by the theory. The maximum cavity width, which 
varies with k, may be taken as a measure of the wake width d'. In 
short, the Riabouchinsky theory, like the notched-hodograph theory, gives 
a flow which depends on the single parameter k and of which the "front 
part" may be used to approximate the flow near a bluff plate. For a 
given base pressure, the drag is very nearly the same as that from the 
notched-hodograph method, as may be expected. The wake width, however, 
is somewhat larger, and the free-streamline shape, of course, is somewhat 
different. The Riabouchinsky theory has been extended to the case of 
wedges by Plesset and Shaffer (ref. 7 ), who also found it necessary to 
use numerical methods to evaluate some of the integrals which occur. It 
could also be adapted to the case of a circular cylinder, as has been 
done here for the notched hodograph. 

In addition, an early paper by Joukowsky (ref. 8) has recently come 
to the author's attention. This gives a general method for the case with 
arbitrary velocity specified along a free streamline . The notched- 
hodograph re suit s for the bluff plate appear there as a special example 
and are interpreted as the flow at a channel mouth which is shielded by 
a flat plate ahead of it . 


WEDGES 


The results of the preceding section give some confidence that the 
notched hodograph will also be suitable for other bluff body shapes and 
that many of the solutions that have been worked out for the Kirchhoff 
condition k = 1 can be generalized in the same way as the case of the 
normal flat plate. The generalization is straightforward for a wedge of 



NACA TN 3168 


15 


arbitrary nose angle 2a , shown in sketch 5* 


For k = 1 the problem 



has been worked out by Bobyleff (see ref. 3, p. 104) . 

The only additional transformation required is one which will open 
the segment in the f^-plane onto the half circle in the £ -plane, from 

where the mapping to the w-plane is identical with that worked out in 
the section "Notched Hodograph." The appropriate transformation is 


£n = 

where 


n = 2a/it 


0 < n 1* 2 


V (15) 


Also, then 
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(The subscripts 1 and n, to distinguish flat plate and wedge, are 
used to prevent confusion in the transformation; they may be omitted 
later.) The mapping from w to £ n then, referring to equation (2), 
is 



The mapping from w to z is 



Since a general solution of the integral in equation (17) could not be 
found for arbitrary values of n, or even for specific values of n other 
than 0 or 1, one case (n = l/2, corresponding to a = 45°) was worked out 
numerically. The resulting values for Cp and d'/d, as functions of k, 
are shown in figure 4 and tabulated in table II. The details for the 
numerical integration are given in appendix B. 


CIRCULAR CYLINDER 


In attempting to apply the free -streamline theory to the circular 
cylinder, two new difficulties are encountered. On the one hand, the 
trace of the cylinder surface in the hodograph plane is not known a 
priori. Second, the separation point in the physical plane is not known, 
as It was In the cases of the plate and wedge. The first difficulty is 
one only of degree — the mapping may always be accomplished in principle. 
The problem of the separation point, however, Is more difficult; it may 
be appreciated from the following discussion. 
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Consider the flow past the curved arc shown in sketch 6. If the 



opening angle 2j3 s is not too large, then separation will occur at S, 

just as in the case of the flat plate. The curvature of the free stream- 

line at separation will depend on the base pressure coefficient Cp g , 

that is, on the parameter k. The lower the base pressure, the more 
sharply will the streamline curve toward the center. As long as the 
"cylinder" consists only of the curved plate shown by the heavy line, 
the stre aml ine can adjust itself to any base pressure, always separating 
at the point S. But if the curved arc Is actually part of a complete 
cylinder, shown dotted, then the radius of curvature of the streamline 
cannot be smaller than that of the cylinder; otherwise the streamline 
would intersect the cylinder. It may, of course, be larger. For one 
particular value of k, that is, of the base pressure, it will be just 
equal to the radius of the cylinder. If the cylinder is actually the 

complete one, then f3 s is net known a priori; but if it be assumed that 

the streamline at separation has the same curvature as the cylinder, then 
there will be a unique value of p s for every value of k. Thus, with 

this assumption, a solution may be obtained, depending as in the previous 
cases only on the parameter k. On the other hand. If the radius of 
curvature at separation Is assumed to be different from (greater than) 
that of the cylinder, then fi 8 will be some other function of k. This 
uncertainty about the conditions at separation makes the free-streamline 
problem of the cylinder considerably more difficult than the case with 
fixed separation points. 


It seems worth while, as an initial step, to work out the case where 
the stre amlin e curvature is the same as that of the cylinder. This had 
already been done for k = 1 by Brodetsky (ref. 9) and later by Schmieden 
(ref. 10). Using a mapping due to Levi Civita, Brodetsky was able to 
obtain an appro ximat e solution by an iteration procedure which converges 
quite rapidly. He found the separation to occur at f3 g = 55 °> giving a 
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drag coefficient Cp = 0*5- It may be expected that for k > 1 the 
values of f3 s and Cjj will be higher. 

To investigate such arbitrary values of k, Brodetsky's method has 
been adapted to the notched hodograph. A more complete discussion of 
the method is given in appendix C. The results are given in figure 5> 
which shows how the separation point (3 S , the drag coefficient C;q, and 
the distance between streamlines d'/d vary with k. Although the 
iteration was carried through only one step, the results up to k = 1.6 
appear to be accurate to a few percent (appendix C). At higher values 
of k the accuracy (with one iteration) becomes more uncertain, so these 
have not been plotted. 

Figure 6 compares an experimental pressure distribution with one 
computed on the basis of the above theory for k = 1.4, chosen to match 
the base pressure on the experimental curve. There is considerable 
improvement over the Kirehhoff case computed by Brodetsky and Schmieden. 
Only in the vicinity of the separation point is there a serious dis- 
crepancy, resulting in quite different values for |3 S in the two cases. 

It is in this region that the uncertainty about the separation condition 
is most noticeable. Clearly the assumption that the streamline has the 
same curvature as the cylinder is ncrt satisfactory. With this assumption, 
there is no adverse pressure gradient (fig. 6), whereas it is well-known 
that to separate a boundary layer on a continuous surface an adverse 
gradient is necessary; it exists in the experimental, case. 

One might, of course, introduce more plausible assumptions about the 
curvature at separation, but for this it appears necessary to go back to 
a study of the boundary- layer separation. That the nature of the boundary 
layer cannot be neglected is clear, since at high Reynolds numbers (above 

105) the separation point moves to p 8 > 90°. The theory in the above 
form is suitable only for values less than ^ 0 ° . 


DISCUSSION 


In each of the examples treated in the preceding sections, and in 
every case to which the theory applies, the solution depends on the single 
parameter k, which cannot be determined without further considerations. 

In this respect the theory is not closed as is the classical (Kirehhoff) 
theory, in which it is simply assumed that k = 1. That choice of k, 
however, is arbitrary, and experience shows.it to be unrealistic. One 
might just as well choose a value in much better agreement with experience. 
For Instance, if the value is determined, experimentally, for one bluff 
body, then it may be used fairly confidently for other shapes at the same 
Reynolds number; at least the result will be" much better than that with 
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k = 1. However, this is a rough observation based on empirical informa- 
tion, not on an essential understanding of the problem. A theoretical 
or semitheoretical basis (e.g., dimensionless analysis) is needed to 
choose k, but the free -streamline theory can take one no further in 
this direction. Some essentially new information must be added; and 
this, it appears, will not be obtained without a consideration of the 
wake, the mechanics of which plays a part in setting the base pressure. 

Due to the fact that the wakes of different bluff bodies develop in 
the same way, from the two separated shear layers, they have many features 
in common. In fact, a wake may be discussed independently of the body if 
its geometrical and velocity scales are known. This is brought out in 
Karmen's analysis of the vortex street, in which the two parameters needed 
to close the problem are a dimension and a velocity (relative to the 
body) . The only function of the body is to determine these two parame- 
ters, or scales; it can hardly have any further influence on subsequent 
developments in the wake (except for Reynolds number effects) . That is, 
the wake is completely determined by a specification of the geometrical 
and velocity scales early in its development, in the transition from the 
body regime to the wake regime. The velocity scale may very appropriately 
be characterized by the velocity along the edges of the free shear layers, 
while the geometrical scale may be specified by the distance between the 
free shear layers. These correspond to the parameters k and d'/d. In 
the free -streamline theory based on the notched hodograph, the relation 
between k and d'/d is determined, for a given body shape, so that 
there Is only one independent parameter, a result which should prove very 
useful . 

While it may be possible to obtain same results, having only the 
geometrical and velocity scales of the wake, it will be necessary even- 
tually to consider the Reynolds number effects. These are related prin- 
cipally to the state of the free shear layers, or of the boundary layer 
before separation. Instead of a Reynolds number based on the body dimen- 
sion, It will probably be more appropriate to introduce the thickness of 
the shear layer and its ratio with respect to d' . This problem has 
received much more attention in the supersonic "base -pressure problem" 
than In the older problem of incompressible flow past bluff bodies. 

Another way In which the results of the notched-hodograph theory 
may be useful Is related to the observations made by Fage and Johansen 
(ref. 2), already mentioned in the "introduction." They observed that 
the frequency of vortex shedding from a bluff body depended not on the 
dimensions of the body but on the distance between the free shear layers, 
and they were able to get a good correlation between bodies of different 
shapes by using this distance in the dimensionless frequency. A single 
parameter like this might be used, in conjunction with the free-streamline 
theory and with measurements of the shedding frequency, to determine the 
drag, for these would give d'/d and thus k and Cp. Here again there 
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is some dependence on Reynolds number, or rather on the ratio of the 
shear-layer thickness to d' . 

Although there seems to be little possibility at present of obtaining 
a theoretical description of the mechanics of the wake, especially in the 
region where it develops from the free shear layers, there is still the 
possibility of finding a correlation between bodies of different shapes 
on a semiempirical basis. Since the completion of this work, a study, 
based on experiment, has been made of the dependence of k on the shear 
layers and on the distance between them. The free-streamline theory is 
combined with some experimental results to obtain a correlation between 
bluff bodies of different shapes, as well as some of the relations between 
wake and body discussed above. The results of this semiempirical study 
are presented in reference 11. 


California Institute of Technology, 

Pasadena, Calif., August 3* 1953- 
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APPENDIX A 
ELLIPTICAL HODOGRAPH 

If the free streamline in the £ -plane is an ellipse (sketch Al), 



Sketch A 1 

the mapping to the w-plane is easily found. (Although the corresponding 
trace of SIS’ in the v-plane (sketch 4 ) is not an ellipse, it will he 
conveniently referred to as the elliptical hodograph.) 

The ellipse in £ is first mapped onto the unit circle, in a new 
plane frcm where the mapping to the w-plane is the same as that in 

the Kirchh'off case of figure 1 . The first mapping is accomplished by a 
Joukowsky transformation. 


i = 


1 + k 
2 




1 - k 
26 i 


while the second is (cf. eq. (2)) 


*1 




These give 
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On the plate 


x = 0 


y = 2^9* + k \j$(l - 9 ) + k tan"- 1 - 


9 


0 < 9 < 1 


while on the free streamline 


x = 1 ^ 9(9 - 1) + logg ^9 - 1 + 


9 £ 1 


y = k | + 2$T 


The breadth of the plate Ib clearly 


d = kit + *(- 


To evaluate the pressure coefficient given by equation (8), 


1= M 0S9S1 


4 


k 2 + ^ ~ ^ 9^1 

9 - 


so that 


s - 


= 1 - 


k 2 9 


(l + k \[l - 9 f 


on the plate and 


C P - 


1 - k c 


1 - k 2 + k^ 
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on the free stre aml ine ■ The average pressure, computed from, equation (10), 
on the front face is 


V = d X, 1 ft " ^ ^ 


= — it + M fc 2 - 1 taDh" 1 — 

die 1 Vk + 1 


which may then be used to find the drag 


Cp = ^pF + k 2 - 1 
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APPENDIX B 
MAPPING THE WEDGE 


The general mapping for the wedge is given by sketch 5 and rela- 
tions (16) and (17). Only the wedge surface and the constant-pressure 
portion of the free stre aml ine will be of particular interest. On these 
w = cp is real and equation (l6) reduces to 



0 < cp < 1 


(Bl) 


£ = e i0 1 ^ <p ^ a^ 


(B2) 


where 


0 = a - ixd 


and 



The average pressure normal to one of the front faces, comparing 
with equation (10), is 




dcp 


(B3) 


where l/q = |M in equation (Bl) . (The subscript has been omitted on 
k n but retained on kq. ) The com ponent in the stream direction contrib- 
utes to the drag the amount 2Cpp sin a. But equation (B3) is averaged 

on L, whereas it is more convenient to compute the drag with reference 
to the base dimensi on d = 2L sin a. On this basis the contribution from 
the front becomes Cpj>, and the drag coefficient is 
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C D = c p F 

1 

Ui 


11 

tfl 

+ - 1 


front face Is 



(*f*A n pi. 

Or ■ 

n—rf 

J 0 ( 

F’ 1+ ! 

)<? " a 2 j 


dcp 


The position of the point B is 


■/; 


Zg = Le^ a + / e^® dcp 


In particular, the distance d' between free streamlines, given by the 
imaginary part of zg, is 


d’ = 2L sin 


a + a 


sin 0 dcp 


■i; 


^a 1 - oa 1- 

= 2L sin a + 2 sin a / cos dp - 2 cos a J sin nui dcp 


and the streamwise distance to B, measured from the base, is 


b = 2 cos 




cos nm dcp + 2 sin 




sin ncu dcp 


The above expressions have been reduced to the operations listed in the 
following summary: 
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g dcp 


H 



1 

g 


dcp 


J 




cos no> dq) 



sin no> dq> 


where 



cd = tan 


-1 
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Then 


L = 



G 



— cot a 
Xj 


b . J 
— = cot a — 
d L 




1 

L 



G 



H 


= 1 


' 2k x 
k-^ 2 + 1 


2n 


H 

G 


C D = + k 2 - 1 



( 


2k, 


, 2n 


ki 2 + 1 , 


H 

G 


The solution for given values of k-]_ and n corresponds to flow over a 
wedge of half -angle a = rut/ 2 and "base -pres sure parameter k = k-J 1 . 
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APPENDIX C 

MAPPING THE CIRCULAR CYLINDER 


Brodetsky (ref. 9) shows how the Kirchhoff flow (k = 1) past a 
curved sire can be computed by means of the mappings shown in sketch Cl. 



Sketch Cl 


First, £ is mapped onto the ft-plane, often used in free -streamline 
theory and defined by 


G = log e 5 = log e i + i0 


(Cl) 


The real part of £1 depends only on the magnitude of the velocity, while 
the imaginary part is the flow direction. The values of 0 on the line CS 
give the flow direction along the surface. For the normal flat plate CS 
in the G-plane would be straight, simply 0 = iit/2 . Now £1 is to be 
mapped onto the w-plane, and the idea used is that the mapping for the arc 
may be obtained by adding a "correction" to the mapping for the flat 
plate, which is known. 
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To accomplish this, 
plane t, defined by 


it is helpful to replace w by the Levi Civita 


1 


£T = 


■» 1 

2iT 


(C2) 


In the r-plane the flow is mapped onto the interior and boundary of the 
semicircle, which has the doublet at I. For the normal plate the mapping 
is 


fl - logg (C3) 

1 - T 

It is assumed that for the curved arc 

ft = log. + fl«( T ) (C4) 

where ft' is to be determined. The singularity occurs in the first term 
so that ft’Cr) may be expressed as a power series, 

A'Ct) = A^t + j A^t^ + . . . (C5) 


The condition that the streamline curvature at separation should be the 
same as that of the surface (see the section "Circular Cylinder") is 
shown by Brodetsky to imply 

A-j- = -1 + aj_ 

A 3 = a l + a 3 
A^ = a^ + a^ 


■^A rather complete discussion of the Levi Civita transformation and 
of the ft-plane, as well as applications to the mapping of curved surfaces, 
is given by Brillouin In reference 12. 
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that is. 


ft = iog e 


1 + T 
1 - T 




(06) 


To evaluate the a^s, the radius of curvature of the surface is expressed 

in terms of them. For a circular arc they are then determined to make the 
radius of curvature constant on the arc. The solution proceeds "by iter- 
ation, starting with a x f 0, a^ = a^ = . . . = 0. Brodetsky showed that 

a-L = 0.0574 gives an arc whose radius is constant within 3^ percent, 
while a-L = O.O585, a^ = -O.OO83 reduces the maximum, discrepancy to 
0.06 percent. The corresponding values of the separation angle P B are 

55. 1° and 55. 0°, while the drag coefficients are 0.4-9 and 0.500, respec- 
tively. This indicates that the convergence is quite rapid, a-j_ already 

giving a fairly accurate re stilt. Therefore, in adapting the method to 
flows where k > 1, it was considered sufficient to use only this one 
term of the iteration. 

For k > 1, the ft k -plane appears as Shown in sketch C2, "notched" 



Sketch C2 


like the corresponding £-plane ; that is, along SB the velocity is con- 
stant (q = 1), while along BI and B'l the streamlines are parallel. Now 
can he mapped onto the 0 -plane (sketch Cl) by the transformation 


sinh^ 


^k “ 


( 


k 2 + 


2k 


4 


sinh c 


a + 



(C7) 


The physical plane corresponding to ftj^ will be called z-^, while that 
corresponding to ; ft is z. 
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Nov the procedure is simply as follows: A value is chosen for a^. 

This determines a certain cylinder shape in the z-plane which corresponds 
to another shape in the zjf-plane. Then k is computed to give a con- 
stant radius of the cylinder in the z k -plane, or as nearly as is possible 
with only the single term a^. Thus for each value of ap there is a 
value of k which gives a circular arc in the z^-plane . The values 

calculated are given in table III. Each arc has a def ini te opening 
angle 20 s which determines the corresponding separation angle 
(sketch Cl). The results, referred to already in the section "Circular 
Cylinder," are plotted in figure 5* Of course, using only a-^ in 

series (C 6 ), it is not possible to make the arc perfectly circular. All 
that can be done is to check the accuracy obtained at each value of k. 

It was found that at k = 1.6 the radius is constant within 2 percent 
and comparable to this at lower values . At higher values of k the 
accuracy becomes lower; therefore the curves of figure 5 have not been 
extended beyond k = 1 . 6 . 

Details of the rather cumbersome calculations for the various quan- 
tities are not shown here. They are analogous to those in Brodetsky’s 
paper, with only the additional introduction of transformation (C 7 ) to 
give the correspondence between the ft- and ft^-p lanes . 

It might be well to point out why it was necessary to introduce this 
transformation from fl^ to ft instead of working directly between ft k 
and t. With the latter procedure the notch would have been distorted, 
since the approximate solution attempts only to obtain a fit on the sur- 
face SC. In the procedure used, on the other hand, the adjustments are 
made in the ft-plane in such a way that the imaginary values on SC are 
constant in the ft^-plane. There is no need then to worry about BIB' 
for equation (C7) maps it exactly onto the notch in the ft k -plane . 


7 
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TABLE I 


FLAT PLATE 


TABLE II 
90° WEDGE 


k 

G 

H 

1.000 

1.106 

1.178 

1.236 

1.360 

1.452 

1.742 

1.733 

1.723 

1.710 

1.696 

0.633 

.639 

.643 

.649 

.660 


9.87 

707 

4.39 


1 

CD 

a* /a 

b/d 

0.637 

.786 

.905 

1.003 

1.261 

1.471 

00 

4 .o 4 o 

2.280 

1.889 

1.438 

1.328 

00 

19.540 

7.800 

3.960 

1.528 

1.036 


TABLE HI 


CIRCULAR CYLINDER 


a l 

k 

0.0574 

1.000 

.0500 

1.175 

,o 4 oo 

1.263 

.0200 

1.387 

0. 

1.470 

-.0400 

1.600 
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Figure 1.- Pressure on flat plate. 
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Figure ij-.- $ 0 ° wedge. 





Figure 5.- -Circular cylinder. 
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